Solitary waves in plasmonic Bragg gratings 
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Light propagation in a Bragg periodic structure containing thin films with metallic nanoparticles 
is studied. Plasmonic resonance frequency, Bragg frequency, and light carrier frequency are assumed 
to be close. Exact solutions describing solitary gap-waves are found, and a light arrest phenomenon 
due to nonlinearity of plasmonic oscillations is studied. 



INTRODUCTION 



The field of photonic crystals, driven by its importance for applications in various areas of photonics, has attracted 
considerable attention during last the ten years [ij- 11 1. The one-dimensional resonant Bragg grating [H-Q or reso- 
nantly absorbing Bragg reflector has been of particular interest. An idealized model of a resonant Bragg grating 
often considers a linear homogeneous dielectric host medium containing an array of thin films with resonant atoms 
or molecules. The thickness of each film is much less than the wavelength of the electromagnetic wave propagating 
through such a structure. The dynamics of ultrashort light pulses in the gratings containing films embedded with 
two- level atoms of the grating has been studied in [l|-0. The transition frequencies of atoms are assumed to be near 
the Bragg frequency. Light interaction with such gratings was described by equations of reduced Maxwell-Bloch type. 
This work demonstrated self-induced transparency in these gratings and the existence of the solitary waves in such 
structures P, |3)l3l- ^^^^ shown Q that bright as well as dark solitons can exist in the spectral gap. 

Progress in nanofabrication of novel optical materials has allowed the design of metal-dielectric nanocomposite 
materials, which have the ability to sustain nonlinear plasmonic oscillations. An example of such a material is a 



dielectric with embedded metallic nanoparticles [13l |- 15|. In this paper we consider ultrashort pulse interaction with 



Bragg structures containing thin films with metallic nanoparticles. The Bragg resonance frequency and carrier wave 
frequency in our considerations are assumed to be close to the plasmonic resonance frequency of the nanoparticles. 
Losses in realistic plasmonic oscillations are of considerable importance. Here we consider an idealized case, wherein 
pulse duration is much shorter than the characteristic time of losses, so that the effects of losses can be neglected. 
In the limit of the slowly-varying envelope approximation we derive governing equations for two counter-propagating 
electromagnetic waves interacting with plasmonic oscillation-induced medium polarization. This system of equations 
represents the two-wave Maxwell-Duffing type model. We find exact solitary wave solutions of this system and, via 
computer simulations, analyze stability of these solutions. 



BASIC EQUATIONS 

We consider a grating consisting of an array of thin films which are embedded in a linear dielectric medium. In 
our derivation of the governing equations we follow [l|-0|, wherein Bragg resonance arises if the distance between 
successive films is a = (A/2) m, m = 1, 2, 3, . . .. 

It was shown in [TtI . [i8| that counter-propagating electric field waves A and B in the slowly-varying envelope 
approximation satisfy the following system of equations: 

<9 1 9\ , q2 d"^ , ^ , 27r(jJo 

where Ago = qo — 27r/a is the mismatch between the carrier wavenumber and the Bragg resonant wavenumber. The 
description of the evolution of material polarization in the slowly- varying amplitude approximation requires modeling 
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of the thin films' response to an external light field. The dielectric properties can be attributed to plasmonic oscil- 
lations, which are modeled b y L orentz oscillators. The simplest generalizations of this model include anharmonicity 
of plasmonic oscillations 3, l3| • In this paper we consider an array of thin films containing metallic nanoparticles 



which haye cubic nonlinear response to external fields [13|, [IJ, [1 



The macroscopic polarization P is goyerned by the equation 

92p dP . ujI 

where ujp is plasma frequency and Ud is the transition frequency between energy leyels resulting from the dimensional 
quantization. Losses of the plasmonic oscillations are taken into account by the parameter Ta- It is assumed that the 
duration of the electromagnetic pulse is small enough that dissipation effects can be neglected. 
Starting from the slowly-yarying envelope approximation, standard manipulation leads to 

+ o)P + |^|P|^7' = -^£„.(.,t). (3) 

Terms varying rapidly in time, which arc proportional to exp(±3zct;ot), are neglected. In this equation £int is the 
electric field interacting with metallic nanoparticles. In the problem under consideration we have Eint — A + B . 

Due to the limitations of nanofabrication, the sizes and shapes of nanoparticles are not uniform. In practice, 
deviation from a perfectly spherical shape has a much larger impact on a nanoparticle's resonance frequency than 
does variation in diameter. This shape variation causes a broadening of the resonance line. The broadened spectrum is 
characterized by a probability density function g{ALu) of deviations Auj from some mean value cUres- When computing 
the total polarization, all resonance frequencies must be taken into account. 

The contributions of the various resonance frequencies are weighted according to the probability density function 
g{Auj); the weighted average is denoted by {V) in equations (HI),©. In what follows, n{ojQ) denotes the refractive 
index of the medium containing the array of thin films, and n„p is the effective density of the resonant nanoparticles 
in films. The effective density is equal to n„p = Nnp{if/a), where iV„p is the bulk density of nanoparticles, if is the 
width of a film, and a is the lattice spacing. 

We study a medium-light interaction in which resonance is the dominant phenomenon. As such, the length of 
the sample is smaller than the characteristic dispersion length. In this case the temporal second derivative terms in 
equations (|f I2p can be omitted. The resulting equations are the two-wave Maxwell-Duffing equations. They can be 
rewritten in dimensionless form using the following rescaling: 

ei ^ A/Ao, 62 ^ B/Ao, p ^ {ATTUJo/[^/eLUpAo])V, ( = {ujp/2c)x, t = t/to. 

Here to — 2y^/ujp, while Aq is a characteristic amplitude of counter-propagating fields. In dimensionless form, the 
two-wave Maxwell-Duffing equations read 

|-^)e2-^e2 = +{ph (4) 

i— + Ap + ^\p\p = -(61-1-62), 

where fj, = {3Ky^/ujoU!p){y/sLUp/ATrLUo)'^AQ is a dimensionless coefficient of anharmonicity, S = 2Ago(c/ti-'p) is the 
dimensionless mismatch coefficient, A = 2y^(LUd — ljq)/liJp is the dimensionless detuning of a nanoparticle's resonance 
frequency from the field's carrier frequency. 

In a coordinate system rotating with angular frequency 5, 

equations (jH) become 

|-|:)./. (5) 

^^ + {A~S)q + ^Aq\'q = -(/1 + /2). 
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Further simplification of the system ([5]) can be achieved by introducing new variables 

/s = -(/l+/2), /a-/l-/2, 

which allow decoupling of one equation from the system of three equations. In these new variables the polarization q 
is coupled with only one field variable. Simple transformations give 

d'fs d^fs „. 5 

z|^ + (A-% + Mkp9 = (8) 
As one can see, we have a coupled system of equations for fs and q. 

SOLITARY WAVE SOLUTIONS 

We consider localized solitary wave solutions of ([51)-® in the limit of narrow spectral line Awg/Aa;^ <C 1, where 
and Acjg are spectral widths of a signal and spectral line g{/S.Lj). In this case the spectral line can be represented 
as a Dirac (5-function: g(Aaj) — 5{/S.lo). Equations ([7]), ^ can then be re-written as follows: 

d^^d^_ dq 

^^ + {A-S)q + ^^\q\\ = (10) 
Scaling analysis of this system shows that solitary wave solutions can be represented as 



2 \ 3/4 



1 / 2v' 
7^ \~v^ 



1/4 

q = qoQn W = ^ ( 7^ ) (v) , (H) 



_ 1 / 2t;2 



77 = (C-vt) 



Here v is velocity of the solitary wave, 77 is a scale-invariant parameter in a coordinate system moving with the solitary 
wave, and functions F^, Qn satisfy the following system of equations: 

-tQ'n + nQn + \Qn\^Qn = Fn. (12) 
The only dimensionless parameter which remains in the system is 



n^iA-6)^'l^, (13) 

which characterizes the deviation of carrier frequency from the plasmonic frequency ujp and the Bragg resonance 
frequency lubt- 

The system of ordinary differential equations (|12p has integral of motion for pulse-like solutions decaying as — > cx3 

W-\Fn\^=0. 
This allows the following parametrization of solutions: 

Fniv) = i?(?7)e^^(''), Qn{v) = R{v)e'^^''\ 
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where R, 0, and ^p are real-valued functions satisfying 

R' = -i?sin(0-i/'), 

0' ^ _cos(0-^), (14) 

ip' + n + R^ ^ cos((/)- 

If we set $ = — -0, then we have 

$'-r2-i?2 = -2cos$, 

i?' = -i?sin$. (15) 

Taking into account equations (|15p wc have the conservation law 

R^ n 

cos$ = h— . (16) 

4 2 

Substituting p6)l into the second equation of ([15]) and subsequent integration gives the following expression for R: 

f7 + 2cosh{/3(ry-?7o)}' ^ ' 



where /3 = \/4 — fi^. The right-hand side is positive real- valued for all 77 if —2 < < 2. Using the conservation 
law (|16p we obtain an expression for <f>: 

$ = 2arctan ( ^ tanh j i/3(7; - ?/o) 1 ) . (18) 



/3 



Now we integrate </>' = — cos $ and find 



17, , /2-n 

- — [1] — i]o) — arctan I — — — tanh 



Finally, we determine tp: 



ip = (f) — 



^ o f2-n , {(3, 

= --(?7 - ??o) - 3 arctan tanh | -(?? - ?/o. 

This pulse exists only if value of the parameter 17 is inside the interval — 2 < 17 < 2. The maximal value of the 
amplitude of this solitary solution is 



A = ^2(2- 17). 

The phases cf) and are nonlinear. Their behavior is asymptotically linear as 77 — > ±00. If 17 = 0, then the limiting 
values of the phases satisfy 

|(/)(oo)-0(-(X3)| = 7r/2, 
IV'(oo) --(/-(-cx))! = 37r/2. 

The total energy of the solitary wave is distributed among the counter-propagating fields and the medium polariza- 
tion. Here we study energy partition among these components. Using equations ([6]), ([7]) and conditions as \ri\ — > cx), 
one can show that 

fa{v)^~-fs{v) = --Fn{v)- (19) 

V V 

We are interested in the energies of the dimensionless fields /i, /2, and of the polarization q. 

h = i(/a-/s) = -y + 

9 = qoQ 
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Finally, for electric fields ei, 62 and for polarization p we have 



2^\l-v^J \v J ^f] + 2cosh{/3(77-77o)} 

e = ^ fJ^y^' (\ i\ (21) 
' 2^V1-^^V V« y V^^ + 2cosli{/3(77-77o)} 



P = 



V7^Vl-«V V^^ + 2cosli{/3(?7-77o)} 



= (5t - -(r/ - - arctan I — tanh j -(ry - 770) 

= f^T - ^(77- 770) - Sarctan ^^-^tanh|^(?7- 770) }> ) . (22) 



The energy of a solitary wave is 



+00 



J |Fj,(77)|2d77 = 8 arctan ^1^. (23) 

— 00 

Finally we have the energies of electric fields and polarization 



= 4 f- - 1 ) (24) 
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Ratios of energies in different fields as well as polarization have the following form 

(25) 
(26) 
(27) 

Therefore, energy partitioning is determined solely by the parameter which is a dimensionless combination of main 
system parameters. 
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NUMERICAL SIMULATION 



The shape and phase of the incident optical pulse are controllable in a real experimental situation. To model 
pulse dynamics in the Bragg grating it is natural to consider asymptotic mixed initial-boundary value problem for 
equations (j4|). We define initial conditions as 

p(C,r)-.0, ei(C,r)^0, e2(C,r)-^0, r ^ -00, (28) 

with no incident field at the right edge of the sample, and with incident field at the left edge defined as follows: 

ei (— 10, r) = ui exp(i6'), 



4.0 exp 



1 /t-3.0 

2 \ 1.5 



= arctan (tanh [1.5 (r — 3.0)]) 



(29) 
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In our case the spatial simulation domain was chosen as [—10,40]. Parameters A — S and fi were given values 

A = 0, (5=1, /i = 1. (30) 

As one can see, we gave the initial pulse the same configuration of phase in a topological sense as would be found in 
a solitary wave solution. This point is important, because otherwise the phase difference cannot relax to the symmetry 
of the stationary wave which is revealed in p^ . As a result, in numerical simulations it is difficult to achieve solitary 
wave type dynamics of the pulse without the proper phase configuration of the incident pulse. 

The second field boundary condition was set to be 

62(40, r)=0. (31) 

In the first experiment we injected a pulse, which is relatively close to the solitary wave solution. The results are 
shown in Figs.[l]l3l which clearly represent two stages of the pulse evolution. In the first stage {t < 7) of evolution we 




FIG. 1: Propagation of pulse. The first experiment. Mapping of the |ei(^,r)| surface. 




FIG. 2: Propagation of pulse. The first experiment. Mapping of the |e2(C,r)| surface. 



observed fast excess energy damping in radiation of quasi-linear waves in both directions and relaxation to a solution 
roughly similar to the stationary solution. Then we observed a stage of pulse shape refinement (7 < t < 30) with 
consequent propagation of the solution very close to (fT7|) . The lower part of the figure [3] shows the creation of multiple 
spatially frozen "hot spots" of the medium polarization. In these hot spots the energy of oscillations is arrested in the 
grating. This "stopping light" phenomenon is due to self-modulation of the medium polarization caused by nonlinear 
effects which we shall discuss in more detail elsewhere. 

During the second experiment we used a pulse of lower amplitude: 



ei(— 10,r) = wexp{i9), 
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FIG. 3: Propagation of pulse. The first experiment. Mapping of the |p(C)''")l surface. 




FIG. 4: Propagation of pulse. The second experiment. Mapping of the |ei(C,r)| surface. 





FIG. 5: Propagation of pulse. The second experiment. Mapping of the |e2(C,''")| surface. 
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FIG. 6; Propagation of pulse. The second experiment. Mapping of the |p((^,r)j surface. 

the second numerical experiment. In this case, computer simulations clearly demonstrate oscillatory behavior of the 
pulse dynamics. Study of the origin of such oscillations will also be presented in a separate publication. In both cases 
we observed a phenomenon similar to self-induced transparency in nanocomposite materials described in [l9| . These 
solitary wave solutions represent the phenomenon of nonlinear light trapping when the wave of medium polarization 
is bound by two optical light fields propagating in the same direction. In the linear limit these two optical fields are 
propagating in opposite directions. 



CONCLUSION 



We studied, both analytically and numerically, solitary waves describing coupled medium polarization and 
light propagation in Bragg gratings with thin films containing metallic nanoparticles. We observed self-induced 
transparency-like phenomena and formation of polarization "hot spots" . The formation of solitary wave solutions 
occurs when energy of an incident pulse exceeds some threshold value. The lower boundary value of this threshold 
can be evaluated as sum of energies (j24p of electric and plasmonic components of coupled wave. 

Stability of nonlinear solitary waves is an important problem which can be addressed by study of collision properties 
and analysis of localized modes of such waves. This problem is the subject of our future investigations. 
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